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This  report  describes  work  in  the  modelling  of  stochastic  phe- 
nomena and  the  development  of  decision-making  techniques  under  risk 
and  uncertainty.  Research  areas  which  received  major  emphasis  were 


(1)  Basic  risk  decision  models,  with  emphasis  on  determining 
the  structure  of  optimal  policies  and  examining  the  impli- 
cations of  different  risk  objectives; 

(2)  Problems  of  data  collection,  estimation,  and  updating  for 


realistic  decision  models. 
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I.  INTRODUCTION 

Almost  all  decision  problems  in  military  and  industrial  operations 
contain  elements  of  risk  or  uncertainty.  These  unknowns  not  only  limit 
the  Information  we  can  gather  about  the  current  state  of  system,  but  make 
future  states  unknowable,  except  in  a statistical  sense,  and  strongly  in- 
fluence attitudes  towards  the  decision  to  be  made.  In  other  words,  the 
utility  of  any  given  decision  can  only  be  forecast  in  a probabilistic  man- 
ner, and  sometimes  cannot  be  evaluated  retrospectively  without  residual 
uncertainty.  Thus,  our  planning  and  control  functions  must  always  reflect 
these  risks  and  uncertainties. 

There  are  a variety  of  models  and  methods  which  appear  under  the 
heading  of  risk  theory;  generally  they  share  the  following  characteristics: 

(1)  There  is  usually  a probabilistic  law  of  motion  governing  the 
system.  In  simple  cases,  this  may  be  a binary  gamble,  increas- 
ing or  decreasing  the  wealth  available  for  future  risk  taking; 
in  more  complicated  models,  a Markov- renewal  or  diffusion  law 
may  be  appropriate. 

(2)  In  addition  to  the  uncertain  dynamics,  there  is  usually  a concern 
about  boundary  conditions,  which  can  lead  to  ruin  (financial,  in- 
surance, gambling  models),  to  catastrophes  (design  of  dams  and 
nuclear  reactors,  responses  to  fires  and  epidemics,  etc.),  or 
simply  to  the  termination  of  the  game  (reliability  and  human- 

or  corporate-lifetime  models).  Sometimes  the  ability  to  delib- 
erately terminate  a (losing)  game  is  the  feature  of  primary  con- 
cern, as  in  optimal  stopping  rules.  Conversely,  in  some  stochastic 
optimization  problems,  there  may  be  some  mathematical  embarrass- 
ment associated  with  continuing  "forever,"  and  turnpike  theorems. 
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discounting,  or  absorbing  states  may  have  to  be  Invoked  to  avoid 
analytic  difficulties. 


(3)  Expected  total  return  may  not  be  a satisfactory  system  objective; 
for  Instance,  the  decision-maker  may  be  just  as  concerned  about 
the  fluctuations  of  reward  under  a certain  policy  as  about  the 
average  reward.  This  may  lead  to  the  use  of  utility  theory  as 

an  axiomatic  way  of  specifying  the  decision-maker's  risk-aversion, 
to  more  explicit  multiple-objective  formulations  or  to  nature-as- 
the-opponent  mlnlmax  functionals.  Many  Interesting  risk-sharing 
problems  require  multiple-player  pareto-optlmallty  approaches. 
Problems  with  long  horizons  usually  require  specific  recognition 
of  the  utility  of  time. 

(4)  Finally,  there  are  always  Important  data  measurement  and  parameter 
estimation  problems  associated  with  decision  under  risk.  Not  only 
may  the  observations  themselves  be  subject  to  error,  but  there  are 
usually  many  more  parameters  to  estimate  for  stochastic  problems 
than  for  deterministic  ones.  An  Important  part  of  many  dynamic 
decision  models  Is  the  provision  for  continuing  updating  of  the 
estimators.  The  mechanics  of  this  are  often  laborious  In  a real 
problem  and  the  costs  of  screening  and  updating  Information  are 
too  often  neglected  In  setting  up  a model. 

During  Che  past  eleven  years  under  this  contract,  "Planning  and  Control 
Under  Risk,"  (DA-31-124-ARO-D-331  and  predecessor  contracts),  we  have  Cackled 
a variety  of  planning  and  control  problems  with  Increasing  emphasis  on  deci- 
sion under  risk;  a list  of  publications  since  1965  Is  In  the  Appendix. 

We  feel  strongly  that  this  will  continue  Co  be  a fruitful  research  area 
In  the  years  to  come.  In  fact,  a heightened  sensitivity  to  the  problems  of 
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risk  within  our  Industries  and  government  agencies  has  been  demonstrated 
In  recent  years,  as  new  problems  of  the  environment,  concerns  over  worker 
and  population  safety,  resource  and  energy  shortages,  etc.  have  been  en- 
countered. Economic  planning  must  now  take  Into  account  possible  sudden 
shifts  In  available  resources,  the  uncertainty  of  long-term  plans,  and  the 
financial  Interdependence  of  world  markets.  Ue  no  longer  have  the  luxury 
of  a certain  world,  but  must  develop  decision  options  for  a variety  of 


alternate  scenarios 
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II.  AREAS  OF  RESEARCH 


A.  Models  of  Decision  Under  Risk 


1.  Optimal  Allocations  In  the  Construction  of  Systems 

In  [1]  a model  of  optimal  allocation  under  risk  was  developed;  the 
problem  Is  to  build  an  n component  system  which  Is  said  to  function  If 
at  least  k (out  of  n)  conq>onents  function.  To  add  a component  we  must 
first  decide  how  much  money  to  allocate  to  that  task;  when  x dollars  are 
Invested  in  a component  then  the  component  will  function  with  probability 
P(x)  . Given  a total  budget  of  A dollars  the  problem  of  Interest  is  to 
determine  how  much  money  should  be  Invested  In  each  component  so  as  to  max- 
imize the  probability  of  attaining  a functioning  system.  The  problem  is 
considered  both  In  the  sequential  and  nonsequential  case,  and  conditions 
under  which  it  is  optimal  to  invest  A/n  dollars  for  each  of  n components 
are  presented.  The  special  case  where  F(x)  min  (x,l)  Is  considered  in 
some  detail.  For  this  case  the  optimal  strategy  is  determined  in  the  se- 
quential case  when  k = 2 and  a conjecture  is  made  in  the  case  of  arbitrary 
k . 

A model  similar  to  the  above  but  with  a different  cost  structure  is 
also  considered  In  [2]  where  once  again  the  problem  is  to  build  an  n com- 
ponent system  which  is  said  to  function  if  at  least  k components  function. 
However  it  is  now  supposed  that  the  components  are  constructed  sequentially 
(with  full  information  as  to  previous  attempts  being  known)  and  if,  after 
constructing  n components,  there  are  only  k - 1 working  ones,  then  a 
penalty  cost  C(l)  is  imposed.  It  is  shown  in  (2)  that  if  C(i)  is  con- 
vex and  if  y^(l)  denotes  the  optimal  amount  to  invest  in  a component  when 
we  have  had  a total  of  k - i successes  in  n - r trials  (the  objective 
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being  to  minimize  the  expected  construction  plus  penalty  costs)  then 

y^(i)  is  nondecreasing  in  i 

and 

y^(i)  is  nonincreasing  in  r 

It  was  then  shown  how  the  above  qualitative  results  could  be  used  to  sub- 
stantially reduce  the  computation  necessary  to  obtain  the  optimal  policy. 

2.  Stochastic  Sequential  Allocation  Models 

In  [3]  a model  in  which  D units  of  resources  is  available  for  in- 
vestment (l.e.,  allocation)  within  a fixed  time  frame  was  considered. 

During  this  time  frame  opportunities  to  Invest  were  assumed  to  occur  in 
accordance  with  a Poisson  process,  and,  as  soon  as  an  opportunity  occurs, 
a decision  as  to  how  much  of  our  available  resources  to  invest  must  be 
immediately  made.  If  y is  Invested  then  an  expected  return  F(y)  is 

gained  and  the  amount  y then  becomes  unavailable  for  future  Investment. 

The  problem  is  to  decide  how  much  to  Invest  at  each  opportunity  so  as  to 
maximize  the  expected  total  return. 

It  was  shown  in  [3]  that  if  F(y)  is  a concave  function  then  if 
y(A,t)  denotes  the  optimal  amount  to  allocate  when  our  present  resources 
is  A , a time  t remains,  and  an  opportunity  is  present  then 

(I)  y(A, t)  is  a nondecreasing  function  of  A , and 

(II)  y(A,t)  is  a nonincreasing  function  of  t . 

These  structural  results  can  be  used  to  simplify  the  necessary  com- 
putation. In  the  special  cases  P(y)  ■ log  y and  P(y)  ■■  y°*  , 0 < a < 1 , 
more  explicit  solutions  for  y(A,t)  were  presented.  When  P(y)  is  con- 
vex it  was  easily  shown  that  y(A,t)  ■ A . 
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3.  Gambling  Models 

In  order  to  obtain  some  Insight  Into  the  structure  of  optimal  policies 
In  risk  models,  a class  of  gambling  models,  useful  as  simple  prototypes  for 
risk  models,  was  considered  In  [4].  For  a variety  of  objectives.  It  was 
shown  that  If  the  game  Is  favorable  to  the  player,  then  he  should  play  as 
timidly  as  possible;  that  Is,  always  make  the  smallest  bet.  A model  In 
which  the  gambler  Is  also  given  the  option  of  working  Is  considered,  and 
It  Is  shown  that  If  the  available  gambles  are  unfavorable  then  the  stra- 
tegy which  minimizes  the  gambler's  expected  time  to  reach  some  preassigned 
goal  Is  the  strategy  that  always  calls  for  working.  For  the  same  model  It 
Is  also  shown  that  If  the  work  option  Is  only  available  at  certain  times 
(namely,  when  the  gambler  Is  broke)  then  the  optimal  gambling  strategy  Is 
to  play  boldly.  These  results  were  obtained  by  developing  some  new  general 
results  In  dynamic  programming,  also  given  In  [4]. 

Currently  under  Investigation  Is  the  problem  where  a decision-maker  Is 
allowed  to  gamble  with  his  objective  being  to  maximize  the  (expected) 
utility  of  his  final  fortune.  It  is  supposed  that  if  he  bets  an  amount  y 
then  his  return  from  this  gamble  will  be  yR  , where  R Is  a random  vari- 
able whose  distribution  Is  known  to  the  gambler,  the  classical  case  being 
when 


(2  with  probability  p 
0 with  probability  1 - p . 

The  objective  of  this  research  Is  to  determine  properties  of  » ^^e 

optimal  bet  when  the  gambler's  fortune  Is  x and  there  are  n gambles  to 
go,  under  the  assumption  of  an  (arlbtrary)  concave  utility  function.  Some 
of  the  Initial  results  obtained  are  that.  In  the  classical  case,  the  opti- 
mal amount  to  save  and  to  strive  for  are  both  Increasing  In  the  gambler's 
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fortune.  That  Is,  In  the  classical  case,  both 


and 

X + yjj(x) 

are  Increasing  functions  of  x . 

4.  Investment  Models 

Many  simple  decision  problems  can  be  paraphrased  In  the  following  way: 
An  entrepreneur,  with  limited  capital.  Is  faced  with  a random  stream  of  va- 
rious Investment  opportunities.  Each  opportunity  Is  a simple  gamble,  that 
Is,  It  requires  a given  Investment  for  a fixed  time  period,  has  a knoim 
probability  of  losing  the  wager,  and  a known  "multiple"  of  gain.  Each 
gamble  Is  usually  superfair.  Although  some  steady-state  optimal,  policies 
are  known  for  this  model,  the  risk  problem  of  Interest  Is  usually  defined 
over  some  uncertain  lifetime,  such  as  the  life  of  an  Individual  or  corpor- 
ation, or  the  time  until  a boundary  (going  broke)  Is  reached.  How  does  the 
desirability  of  the  investment  vary  with  size  of  capital,  the  average  re- 
maining horizon,  and  investment  stream  fluctuations? 

Another  Investment  model  of  Interest  is  a consumption-savings  scheme 
In  which,  at  the  beginning  of  each  time  period,  a decision  Is  made  as  to 
how  much  of  one's  available  capital  is  to  be  consumed  (the  remainder  to  be 
Invested).  The  Interest  rate  for  savings  changes  stochastically  from  period 
to  period.  While  previous  models  (see  [6]  for  relevant  references)  seem  to 
suppose  that  the  decision-maker  does  not  know  the  Interest  rate  at  which 
his  savings  will  grow  when  the  decision  Is  made  (the  Interest  rate  Is  as- 
sumed to  be  a random  variable  whose  value  Is  found  out  after  the  decision) 
we  propose  to  study  the  case  where  the  value  of  the  Interest  rate  Is  known 
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before  the  decision  is  made.  Thus,  for  example  our  state  space  now  consists 
of  a three-dimensional  vector  (n,x,r)  where  n is  the  number  of  time 
periods  remaining  In  the  problem,  x Is  the  present  fortune,  and  r Is 
the  present  Interest  rate.  Another  model,  with  possibilities  for  buying 
Insurance  against  risk,  is  described  in  the  next  sector.  Of  course,  models 
in  the  investment  area  are  closely  related  to  the  gambling  models  described 
In  the  last  section. 

5.  Risk-Sharing  Contracts 

In  many  situations,  it  is  possible  to  reduce  a risk  by  paying  a premium 
to  a third  party,  such  as  an  insurance  company,  or  by  trading  similar,  but 
more  nearly  uncorrelated  risks  with  other  risk- takers.  If  such  an  arrange- 
ment continues  over  many  years,  there  are  a variety  of  "fair"  arrangements, 
taking  into  account  the  differing  utility  functions  of  the  parties  involved, 
the  differing  risks,  the  need  for  capital  reserves  (to  avoid  ruin),  etc. 

Thus  a variety  of  secondary  criteria  must  be  invoked  to  decide  between  these 
different  "fair"  arrangements.  Most  of  the  current  literature  is  in  the 
actuarial  journals,  with  particular  emphasis  on  the  problems  of  reinsurance 
corporations.  However,  the  principles  involved  are  applicable  to  any  busl-  ; 

ness  arrangement  where  liability  for  financial  risks  may  be  shared. 

A recent  thesis  [6]  explores  in  detail  the  trade-off  between  payout 

1 

(in  the  form  of  dividends)  and  hedging  (in  the  form  of  insurance)  which  is 

i 

available  to  a corporation  (here  described  as  an  insurance  company) . In  j 

one  part  a dynamic  programming  model  is  formulated  for  the  dynamic  dividend 

— (re) insurance  decision  problem,  and  closed-form  optimal  policies  are  found 

for  a particularly  interesting  class  of  utility  functions;  this  model  is  j 

similar  to  those  described  in  the  last  two  sections.  The  other  part  of  the  I 


thesis  considers  what  happens  when  several  such  corporations  meet  to  arrange 
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a mutual  pooling  of  risk.  The  existence  of  an  equilibrium  price  function 
in  the  resulting  (re) Insurance  market  is  then  demonstrated. 

Of  current  Interest  are  problems  in  which  risk  contracts  can  be  broken 
at  any  time,  and  neither  party  has  complete  knowledge;  information  updating 
occurs  via  linear  Bayesian  methods  (see  next  section  on  credibility  theory). 
The  definition  of  what  constitutes  a "fair"  dynamic  risk  premium  is  of  cen- 
tral importance  here. 

B.  Estimation  Problems 

Bayesian  Estimation 

Most  decision  problems  under  uncertainty  require  processing  of  large 
amounts  of  data,  either  for  initial  estimates  of  the  parameters,  or  to  pro- 
vide continuing  updating  as  new  information  is  received  from  prior  decisions. 
Thus,  efficient  methods  of  estimation  are  Important. 

Suppose  we  can  observe  a random  variable  x which  depends  upon  a 
parameter  0 in  a known  way,  via  the  tvketyhood  density  p(x|0)  ; we  as- 
sume a prior  distribution,  u(6)  , is  available.  Prior  predictions  about 
an  average  x can  then  be  made  through  the  mixed  density  p(x)  ■ E„p(x|6)  . 
Now  suppose  that,  as  a result  of  some  decision  about  experimentation  and 
sampling,  we  can  observe  n Independent  samples  of  the  random  variable, 

X ■■  {x^  ••  x^  ; t ••  1,2,  ...,  n}  . By  using  Bayes'  theorem,  the  density  of 
6 is  updated  as  follows: 

np(Xj.|e)u(0) 

" /np(Xj.U)u(.p)d<|/ 

In  decision  problems,  we  are  not  usually  so  Interested  in  forecasting  the 
parameter  0 , as  we  are  in  the  foreoaet  density  for  S' 


y , the  next 
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observation.  Thus,  practical  control  problems  require  the  calculation  of: 

(2)  p(y|x)  “ » J*p(y|e)u(e|x)de 

or  Its  moments.  For  Instance,  In  decision  analysis,  as  decisions  to  perform 
experiments  are  made,  the  Information  from  the  experiments  Is  used  to  recom- 
pute expected  utilities  of  one  or  more  future  actions. 

Practically  speaking,  computations  via  (1)  and  (2)  are  laborious,  and 
require  either  large  computer  capability,  or  the  use  of  a few  well-known 
natural  conjugate  prior  families  of  llkelyhood  and  prior.  This  requires 
knowing  (or  being  able  to  assiune)  a great  deal  of  structure  Information 
about  u(6)  and  p(x|9)  . But,  as  Indicated  above,  most  of  this  Informa- 
tion Is  "wasted,"  especially  If  we  only  want  moments  or  an  expected  utility. 

Credibility  Estimation 

In  the  1920* s American  actuaries  developed  a class  of  estimators  called 
"credibility  formulae"  to  predict  the  mean  of  the  next  risk  outcome,  > 

given  the  "experience  data,"  x . Using  simple  heuristic  arguments,  they 
derived  a linear  forecast  formula  f(x)  : 

E{x^llx  }S:f(x)  = (1  - z)*m  + I x^ 

(3) 

n 

^ “ n + N 

1 V 

where  m Is  the  prior  mean  of  p(y)  (no  data),  L the  sample  mean 

of  data,  and  N Is  a time  constant,  chosen  originally  In  an  ad  hoc  manner. 
The  most  Interesting  features  of  this  formula  are:  (1)  It  uses  the  data  In 

a simple  linear  fashion,  mixing  It  with  the  prior  estimate;  (2)  It  worked 
extremely  well  for  40  years  In  casualty  Insurance  experience  rating. 
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In  the  1950' 8 It  was  discovered  that  (3)  was  in  fact  the  exact  Bayesian 

resxilt  E(x  lx)  (gotten  from  (1)  and  (2))  If  the  prior  and  likelyhood 
n+l  — 

were  the  natural  conjugate  priors:  Beta-Blnomlal,  Gamma-Folsson,  and  Nor- 

mal-Normal. Then  In  1967,  It  was  shown  that  the  credibility  forecast  (3) 
is  the  best  least- sqwxrea  approximation  for  arbitrary  distributions,  if 
N is  chosen  as  the  ratio  of  the  two  components  of  prior  variance. 

Recent  Research 

In  a research  effort  supported  by  ARO  since  1973,  we  have  been  able  to 
find  a variety  of  interesting  and  useful  extensions  to  credibility  estima- 
tions. In  [10]  [12],  the  class  of  llkelyhoods  for  which  (3)  is  exact  was 
extended  to  all  distributions  where  the  range  of  x does  not  depend  upon 
6 ^ and  for  which  the  sample  mean  is  (one  of)  the  sufficient  statistics. 

It  is  well  known  that  this  implies  p(x|0)  belongs  to  the  Koopman-Pltman- 
Darmois  exponential  family ^ which  includes  a very  large  number  of  statis- 
tically important  distributions,  and  for  which  the  natural  conjugate  prior 
families  is  easy  to  obtain.  The  forecast  (3)  only  requires  estimation  of 
the  prior  mean  and  two  components  of  variance  in  x-space;  the  prior  itself 
does  not  enter  explicitly.  Thus,  credibility  forecasts  have  very  valuable 
properties  for  forecasting  average  values: 

(1)  They  are  almost  "distribution-free,"  requiring  the  estimation  of 
very  few  prior  moments,  always  in  the  space  of  observables; 

(2)  The  data  are  used  in  a simple  way,  as  a sample  mean  or  other  suf- 
ficient statistic. 

Furthermore,  the  credibility  approach  can  be  generalized.  If  other 
sufficient  statistics  are  known  (such  as  the  2nd  moment)  to  be  appropriate, 
it  is  pointed  out  in  [9]  that  a formula  similar  to  (3)  also  applies  to  these 
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generalized  moments. 

The  basic  model  has  also  been  extended  to  the  case  where  x is  a veo- 
tor-vatued  random  variabley  whose  prior  distribution  is  to  be  updated  using 
a matrix  X of  observed  data.  Reference  [7]  gives  the  least- sqxiares 
linearized  approximation  for  the  Bayesian  forecast;  reference  [10]  shows 
that  the  forecast  is  exact  for  multi-dimensional  exponential  families,  such 
as  the  multinomial  and  multinormal  distributions.  The  resulting  forecast 
formula  is  a vector  version  of  (3) , requiring  only  prior  estimates  of  the 
(vector)  mean  and  (matrix)  covariances. 

In  another  direction  [8],  it  has  been  shown  that  a credibility  formula 
can  be  obtained  for  an  estimate  of  the  probability  P(ylx)  “ PrCx^j^  ± y|x] 
for  a fixed  value  of  y . Updating  Involves  mixing  one's  prior  estimate  of 
this  probability  with  the  sample  frequency  of  success  in  a form  similar  to 
(3),  instead  of  direct  conq>utatlon  via  (1)  and  (2).  This  procedure  is  ex- 
act in  only  one  trivial  case  (Beta-Bemoulll  distributions) , but  computer 
simulations  have  shown  that  it  is  a good  approximation  to  the  exact  Bayesian 
result,  being  unbiased  and  having  a variance  less  than  twice  the  theoretical 
minimum.  And  its  simplicity  recommends  itself  for  updating  models  with  a 
large  data  base. 

Paper  [13]  explores  a large  class  of  model  extensions  which  are  of 
practical  use,  such  as  special  data  structures,  correlated  observations, 
evolutionary  models,  and  integration  of  cohort  data.  When  this  paper  ap- 
peared, the  Principal  Investigator  was  host  to  an  Actuarial  Research  Con- 
ference on  Credibility  Theory  at  Berkeley,  sponsored  by  the  American  Society 
of  Actuaries  and  Casualty  Actuarial  Society.  Interest  in  credibility  theory 
thus  continues  high  within  actuarial  circles,  in  addition  to  the  many  dif- 
ferent applications  possible  in  industrial  and  military  estimation  problems. 
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The  next  five  papers  [14]  [15]  [16]  [17]  [18]  appeared  while  one  of 
the  Investigators  was  on  sabbatical  leave  at  the  International  Institute 
of  Applied  Systems  Analysis  at  Laxenburg,  Austria,  from  September  1974  to 
June  1975.  [14]  analyzes  the  two  classes  of  covariance  matrices  for  which 

one  Is  able  to  get  closed-form  solutions  to  least-squares  problems.  [15] 
continues  the  work  of  [13]  on  collateral  data,  and  develops  a hierarchical 
model  which  explains  how  related  data  structures  may  be  Integrated  In  a 
single  forecast.  [16]  develops  Bayesian  credibility  estimators  for  simple 
Inverse  regression  problems  which  arise,  for  example.  In  the  calibration  of 
Instruments  and  other  simple  linear  laws.  [17]  provides  a unified  overview 
of  credibility  as  applied  to  general  Bayesian  regression  problems,  based 
upon  Initial  work  by  Hachmelster  «md  Taylor  which  was  given  at  the  Actuarial 
Research  Conference  on  Credibility;  this  approach  has  Important  ties  with 
the  Kalman  filtering  methods  used  In  statistical  control  theory.  [18]  ex- 
plores the  application  of  credibility  to  a problem  In  material  verification, 
such  as  might  arise  In  a Safeguards  program  for  radioactive  material  account- 
ability. Because  of  the  relative  Inaccessibility  of  these  papers  to  Army 
personnel.  It  Is  planned  to  reissue  them  In  O.R.  Center  format  In  the  near 
future,  as  soon  as  proper  permission  can  be  obtained. 

An  Important  problem  In  any  forecast  Is  the  estimation  of  trend  com- 
ponents. Although  this  can  be  easily  formalized  In  any  Bayesian  regression 
problem  by  expanding  the  number  of  unknown  coefficients.  [19]  Indicates  how 
credibility  trend  analysis  can  be  carried  out  In  a compact  and  efficient 
manner  by  using  the  structure  explicitly  to  reduce  the  problem  dimensions. 

Paper  [20]  explores  the  application  of  credibility  theory  to  the  Bayesian 


estimation  of  network  flows,  which  arise  In  various  logistic  and  material  ac- 
countability problems.  This  problem  la  notoriously  "unidentifiable"  within 
the  framework  of  classical  estimation  theory,  and  a full  Bayesian  analysis 
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vrould  require  very  large-dlmenslonal  multi-dlmenslonal  priors  and  likely- 
hoods  to  be  specified. 

Also  part  of  this  research  effort  is  [21],  ^ich  develops  an  exact 
Bayesian  model  for  life-testing  models.  In  which  many  of  the  observations 
are  of  the  {SP  > x}  type.  This  report  was  supported  by  ARO  under  Contract 
DAA629-76-6-0042  with  Richard  Barlow. 
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